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O ; ABSTRACT 

\^ • The dynamical ages of the opposite lobes of selected giant radio sources are estimated 

' using the DYNAGE algorithm of Machalski et al., and compared with their spectral 

Q^, ages estimated and studied by Jamrozy et al. in Paper II. As expected, the DYNAGE 

fits give slightly different dynamical ages and other model's parameters for the oppo- 
site lobes modelled independently each other, e.g. the age ratios are found between 
~1.1 to ~1.4. Demanding similar values of the jet power and the radio core density 
for the same source, we look for a self-consistent solution for the opposite lobes, which 
results in different density profiles along them found by the fit. We also show that a 
departure from the equipartition conditions assumed in the model, justified by X-ray 
observations of the lobes of some nearby radio galaxies, and a relevant variation of the 
magnetic-field strengths may provide an equalisation of the lobes' ages. A comparison 
of the dynamical and spectral ages shows that a ratio of the dynamical age to the 
spectral age of the lobes of investigated giant radio galaxies is between ~1 and ~5, i.e. 
, is similar to that found for smaller radio galaxies (e.g. Parma et al. 1999). Supplement- 

• ' ing possible causes for this effect already discussed in the literature, like uncertainty of 

assumed parameters of the model, an influence of a possible departure from the energy 
equipartition assumption, etc., the further two are pointed out and discussed: (i) a 
difference between the injection spectral indices describing the initial energy distribu- 
tions of the emitting relativistic particles determined using the DYNAGE algorithm in 
the dynamical analysis and in the classical spectral-ageing analysis, and (ii) a different 
influence of the axial ratio of the lobes in estimation of the dynamical age and the 
spectral (synchrotron) age. Arguments are given to suggest that DYNAGE can better 
' take account of radiative effects at lower frequencies than the spectral-ageing analysis. 

The DYNAGE algorithm is especially effective for sources at high redshifts, for which 
an intrinsic spectral curvature is shifted to low frequencies. 
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1 INTRODUCTION 

There are several approaches to estimate the age of a classi- 
cal double radio source, beginning with a use of the projected 
linear size of that source and estimating the speed of its ex- 
pansion to that size. Direct measurements of the expansion 
speed inferred from the proper motions of the hot spots in 
compact symmetric radio sources gave values of about 0.2c 
- 0.3c (cf. Owsianik, Conway & Polatidis 1998; Owsianik 
& Conway 1998). An extrapolation of these motions back 
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in time indicates very young ages of such sources, being of 
the order of 10^ - 10* years. On the other hand, all the 
analytical models of the dynamics and radio-emission prop- 
erties of powerful double-lobed radio sources (e.g. Scheuer 
1974; Begelman & Cioffi 1989; Falle 1991; Nath 1995; Kaiser, 
Dennett-Thorpe & Alexander 1997; Blundell, Rawlings & 
Willott 1999; Manolakou & Kirk 2002; Kino & Kawakatu 
2005) predict that those speeds reduce gradually by an or- 
der of magnitude or even more with the source age. 

The above has been confirmed by the classical spectral- 
ageing analysis. There is no doubt that radio continuum 
spectra in different parts of an extended radio source con- 
tain important information about the various energy losses 
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and gains of the radiating particles during the hfetime of 
the source. With the assumptions that (i) these particles are 
immersed in a uniform magnetic field, (ii) they are not sig- 
nificantly reaccelerated within the source (their lobes), and 
(iii) there arc no significant mixing of new and old particles 

- the observed radio spectrum should steepen with increas- 
ing distance from the place of the last acceleration, i.e. from 
the hot spots. These predictions have been detected in many 
radio sources and used to estimate the radiative ages of the 
emitting particles and the expansion speeds in several sam- 
ples of powerful SCR sources (e.g. Myers & Spangler 1985; 
Alexander & Leahy 1987; Leahy, Muxlow & Stephens 1989; 
Carilli ct al. 1991: Liu, Poolcy & Riley 1992), in samples of 
low-luminosity and medium-luminosity radio galaxies (e.g. 
Klein et al. 1995; Parma et al. 1999), as well as in samples 
and/or of individual 'giant'-sized radio sources (e.g. Lacy et 
al. 1993; Saripalh et al. 1994; Mack et al. 1998; Schoenmak- 
ers et al. 1998, 2000; Lara et al. 2000). 

However, the observed steepening of the spectrum need 
not be entirely due to radiative energy losses. A possible 
evolution of the local magnetic fields, a bulk backflow and 
significant mixing of the lobe material, or the difficulties in 
disentangling the effects of the various loss processes have 
been pointed out in a number of papers (cf. Rudnick, Katz- 
Stone & Anderson 1994; Eilek & Arendt 1996; Jones, Ryu 
& Engel 1999). The spectral steepening due to the energy 
losses is parameterised by a single 'break' in the observed 
spectrum, !/br. The spectral age of the particle population 
with the observed Z4,r within a constant magnetic field of 
strenght B is proportional to B~^^'^v^^^^ . If B decreases 
while the source (its lobes) expands, the spectral age should 
overestimate the true age of the source. However this age is 
usually found to be lower than the source's age inferred from 
the dynamical considerations (cf. Kaiser 2000, hereafter re- 
ferred to as K2000). 

For the first time the problem of how to reconcile the 
spectral and dynamical ages was undertaken by Blundell & 
Rawlings (2000). They discussed how these two can be quite 
discrepant from one another rendering use of the classical 
spectral ageing method inappropriate. Moving beyond the 
traditional bulk backflow picture and considering alternative 
means of the transport of high-energy particles, the authors 
explained the spectral steepening along the lobes not pre- 
dominantly by synchrotron ageing but by gentle gradients 
in the magnetic field. They contended that spectral ages can 
give meaningful estimates of dynamical ages only when these 
ages axe less or much less than lO'^ years. The same prob- 
lem was studied by K2000 who extended the spectral-ageing 
methods including the underlying source dynamics into the 
age estimates. The author claimed that if the bulk backfiow 
and energy losses of the relativistic electrons, both radia- 
tive and adiabatic, are self-consistently taken into account, 
the discrepancies between spectral ages and dynamical ages 
arising from the earlier methods can be resolved. However 
analysing the K2000 model, Machalski et al. (2007) realised 
that in a majority of the extended FRII-type radio sources 

- even those without distorted lobe structures - the surface- 
brightness profiles are far from the expected smooth shapes, 
making the fitted free parameters of the model highly un- 
certain. Besides, the K2000 method requires rather high- 
resolution observations of the radio lobes. But usually high- 
resolution observations of low-brightness sources cause a se- 



rious loss of the fiux density, so that the K2000 method 
can in fact only be applied to the strongest sources such as 
CygA. 

In Jamrozy et al. (2008, Paper II of this series), multifre- 
quency observations with the Very Large Array (VLA) and 
the Giant Metrewave Radio Telescope (GMRT) (Konar et 
al., 2008, Paper I of this series) have been used to determine 
the spectral ages of ten selected giant-sized radio galaxies 
(GRGs). Using the classical spectral-ageing approach and 
applying two different formulae for the cquipartition mag- 
netic field estimates: the classical formula of Miley (1980) 
and the revised formula of Beck & Krause (2005) - the spec- 
tral age distributions along the main axis of the lobes were 
analysed. In that paper we found that practically the statis- 
tics of the derived age is independent of the equipartition 
field formula applied, though those ages can be quite dif- 
ferent for individual sources. We also analysed the injection 
spectral indices characterising an initial power-law energy 
distribution of the emitting particles. Those indices, deter- 
mined by fit to the observed radio spectra with the SYNAGE 
algorithm of Murgia (1996) were found to be correlated ei- 
ther with luminosity or redshift, as well as with linear size 
of the source. 

In this paper dynamical ages of the lobes of the ten 
GRGs studied in Paper II are determined and compared 
with their spectral ages. Using the DYNAGE algorithm 
(Machalski et al. 2007), we derive the age, the average ex- 
pansion velocity of the lobes' heads, the effective injection 
spectral index which approximates the initial electron con- 
tinuum averaged over a very broad energy range and over 
the present age of source, and other dynamical properties 
of the sources like their jets' power, central density near the 
radio core which determines the local environment density 
in which the jets propagate, and the internal pressure in 
the lobes. In Section 2 we describe briefly the DYNAGE al- 
gorithm for fitting the dynamical parameters of the model 
to the observational data given in Section 3. In Section 4 
we present the results of the fit, while the discussion of the 
results, especially of the two factors causing a difference be- 
tween the dynamical and the spectral ages are given in Sec- 
tion 5. 



2 THE DYNAGE ALGORITHM 

This algorithm is an extension of the analytical model for the 
evolution of FRII-type radio sources combining the dynam- 
ical model of Kaiser & Alexander (1997) with the model for 
expected radio emission from a source (its lobes or cocoon) 
under the influence of energy loss processes published by 
Kaiser, Dennett-Thorpe & Alexander (1997, hereafter KDA 
model). One of the basic assumptions of the KDA model 
is a continuous delivery of kinetic energy from the active 
galactic nucleus (AGN) to the radio lobes through the jets 
(e.g. Falle 1991). The jets terminate in strong shocks where 
the jet particles are accelerated and finally infiate the co- 
coon. The density distribution of unperturbed external gas 
is approximated as 

p(r) = po (—] ; for r > ao, (1) 
vao/ 
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where po is the central density at the core radius ao, and 
the exponent /3 describes the density profile in the simpli- 
fied King's (1972) model. Such a distribution of the ambi- 
ent medium is Eissumed to be invariant with redshift. As in 
KDA, the lobe (i.e. half of the cocoon) is approximated by a 
cylinder of length D and base diameter b, so that its volume 
Vc is determined via the axial ratio, Rt = D/h. The cocoon 
expands along the jet axis driven by the hotspot plasma 
pressure Ph and in the perpendicular direction by the co- 
coon pressure pc- The ratio of these pressures is constant for 
a given radio lobe and depends on its Rt. In our calculations 
we use the empirical formula taken from K2000 

Phc = Ph/Pc = (2.14 - 0.52/3)i4°^-°-^5''. 

Thus the model predicts a self-similar expansion of the co- 
coon (lobe) and gives analytical formulae for the time evo- 
lution of its geometrical and physical parameters, e.g. the 
length of the lobe 



1/(5-/3) 

L{t) = ci ( ^ ) t"^'-'\ 
and the cocoon pressure 



(2) 



(2-/3) 
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where ci is a dimensionless constant, Fx - the adiabatic in- 
dex of unshockcd medium surrounding the cocoon (lobe), t 
- the time elapsed since the jet started from the AGN, i.e. 
it is the source's (here its lobes') actual age. The cocoon's 
pressure determines the energy density within it via the adi- 
abatic index Fc of the cocoon as a whole, depending on the 
relative pressures of relativistic electrons, thermal particles, 
and magnetic 'fluid' 

«c(t) =Pc(t)/(rc~i), 

where Uc = -ub + We ( 1 + fe' ) , with the equipartition condition 

WB 1 -h P Oinj + 1 



(4) 



where fc' is the ratio of the energy density of thermal par- 
ticles to that of the electrons, and p is the power index of 
the power-law energy spectrum of radiating particles. The 
magnetic field (assumed to be completely tangled) with the 
energy density ub and adiabatic index Fb satisfied the rela^ 
tion 



1/9 r (4+/3)rB ] 



(5) 



As in KDA, we assume that the jets consist of an electron- 
positron plasma and that the energy contribution from rel- 
ativistic protons is completely negligible. 

The above dynamical equations are supplemented with 
an integral giving the total radio emission from the cocoon 
at a firequency v. The cocoon is split into many small volume 
elements 5V, each of which is allowed to evolve by expand- 
ing adiabatically, i.e. changing the pressure from the hotspot 
pressure Ph(ti) to the cocoon pressure Pc(ti) as a function of 
injection time U. Tracing the effects of adiabatic expansion, 
synchrotron losses (with the assumed effective isotropisa- 
tion of the particles' pitch angle distribution), and inverse 



Compton scattering on the cosmic background radiation in 
the volume elements independently, the radio power of the 
cocoon at a fixed observing frequency is obtained by sum- 
ming up the contributions from all elements, resulting in the 
integral over U (equation (16) in KDA). The integral is not 
analytically solvable and has to be calculated numerically. 

The DYNAGE algorithm allows one to determine the 
values of four of the model's free parameters, i.e. the jet 
power, Qjet, central core density, po, injected spectral index, 
Oinj, and dynamical age, t. The determination of their values 
is possible by the fit to the observational parameters of a 
source: its projected linear size, D, the volume of the cocoon, 
Vc, the radio luminosity, P^, and the radio spectrum, a^, 
which provides P^j at a number of observing frequencies 
i = l, 2, 3.... 

The values of several free parameters of the model have 
to be assumed. The assumed values are listed in Table 1, 
where 7i,min and 71, max arc the Lorentz factors determining 
the energy range of the relativistic particles used in integra- 
tion of their initial power-law distribution. Since there is a 
mixture of three different fluids in the cocoon, this is likely 
that its equation of state varies within the cocoon and with 
time. In the present calculations we choose the KDA 'Case 
2' where the cocoon's mixed material has a non-relativistic 
equation of state but the energy density of the magnetic field 
is proportional to that of the relativistic particles. As dis- 
cussed in KDA, in this case the value of ^ increases with time 
which implies that after some time ub starts to dominate the 
total energy density in a given volume element 5V . However 
at the same time the total volume of the cocoon increases 
by a much larger factor, so that the approximation of the 
adiabatic index of the whole cocoon of 5/3 is justified. The 
calculations show that the model fits obtained either with 
Fc=5/3 or Fc=4/3 is negligible if k' is kept constant. More 
pronounced is influence of a value of k' itself. An increase of 
fc' from zero to ten causes a ~10% to ~20% increase of the 
age due to a thermal expansion of the lobes (for the largest 
values of Rt it can be more than 30%, cf. Machalski et al. 
2008). This effect, though ruling the entire time scale, does 
not affect age differences between the opposite lobes (as far 
as fc' values are comparable in both lobes). However, with- 
out an independent evidence about the fractional content 
of thermal particles, holding its value as zero seems to be 
acceptable when we analyse the above differences. 

Furthermore, we assume the orientation of the jet axis 
to the observer's line of sight, 9. Because of the giant lin- 
ear sizes of the investigated sources, we assume 9 = 90° 
for almost all the lobes except two for which the observed 
asymmetries in the lobes' separation and/or in their bright- 
ness suggest 9 < 90°. The details are given in Section 3. An 
extensive discussion of limitations of the DYNAGE method 
and dependence of the age solution on the assumed values 
of the model's free parameters are given in Machalski et al. 
(2007). For example, changing the value of 7i.min from 1 to 
10 decreases the age of a source (lobe) by about 2 %-5 %, 
while varying the value of ao between 20 kpc and 5 kpc 
changes the age from about —10% to about -|-20%, which 
is not larger than the age uncertainties determined for the 
GRGs analysed in this paper (cf. Table 3). 
The fitting procedure consists of three steps: 

(1) For a given value of Qinj and a number of values of 
t, the values of Qjet{cniij,t, P^) and po(a!inj, Pi^) are deter- 
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Table 1. Assumed values of the model free parameters 



ao /3 7i,min 7i,max 1-cilx l-B ^ " 

lOkpc 1.5 1 lO'^ 5/3 4/3 90° 



mined by the fit of equation (2) to the deprojected linear size 
of a lobe D/ sin 9 and equation (16) of KDA to the lobe's 
luminosity at a given frequency P^. 

(2) Performing the above procedure for all of the lumi- 
nosities chosen to represent the lobe's spectrum, we get a set 
of solutions for Qjet and po shown in Fig. la. As the source 
(its lobe) must have the same age at any of the observing 
frequencies, we search for an age at which the parameters 
Qjet and po have possibly identical values. A perfect inter- 
section of all the Qjet — Po curves corresponding to differ- 
ent observing frequencies at some particular age is expected 
if the observed spectrum agrees with the theoretically pre- 
dicted one in the framework of the 'continuum injection' 
(C.I.) model of energy losses (cf. Myers & Spangler 1985; 
Carilli et al. 1991). However, the observed spectrum of the 
analysed sources (and their lobes) may depart from such a 
theoretical shape, therefore the 'goodness' of the intersection 
is quantified by the A{t) measure defined by equation (4) in 
Machalski et al. (2007). A minimum of A(t) (i.e. 'the best' 
intersection of the Qjct — po curves for different frequencies 
and for the given value of ainj ) is considered as an estimate 
of the source's (lobe's) dynamical age t(ainj) (Fig. lb). This 
minimum also distinguishes a value of Qjet(Q:inj) (and po, Pc, 
Wc, all dependent on a value of Oinj). 

(3) The steps (1) and (2) are repeated for a number 
of Oiinj values. By varying this parameter during the fitting 
procedure, we search for its value which provides a minimum 
of the product Qjet x t, i.e. a minimum of the kinetic energy 
delivered to the given lobe during the age t. Following the 
original KDA assumption during integration of the expected 
radio emission, the minimum energy condition is initially 
fulfilled in each volume element of the cocoon, therefore the 
minimum of Qjot x t corresponds to a minimum of the total 
energy density in it. 



3 THE DATA 

The observational parameters of the GRGs under investiga- 
tion and their lobes are given in Table 2. Columns (1), (2), 
(3), and (7) are self explanatory. Columns (4) and (8) give 
linear size D of the opposite lobes (the first line) and their 
axial ratio Rt (the second lino); both values with their stan- 
dard error. The value of Rt is determined on the maps with 
the highest sensitivity to surface brightness available as the 
ratio of D and the largest deconvolved width of transver- 
sal cross-sections through the lobe, b. The standard errors 
of D and _Rt result from the uncertainties of angular size 
of the lobe's length and its base diameter. Columns (5) and 
(9) give the observing frequencies selected from the available 
spectral data and used to calculate the radio luminosities, 
P,/, necessary for the fitting procedure (cf. Section 2). For a 
technical reason, we limited the fitting procedure to four fre- 
quencies only, possibly spanning the largest frequency range. 
The letter preceding the observing frequency indicates the 
radio telescope (array) or the radio survey with (or in) which 



a relevant flux density was measured: G - for GMRT, V - 
for VLA, and W - for WSRT. Mostly the flux densities pub- 
lished in Paper I are used, but when necessary, the spectral 
data are supplemented with flux densities taken from other 
observations (surveys): VLSS (Cohen et al. 2007), 7C (Poo- 
ley et al. 1998; Riley et al. 1999), WENSS (Rengelink et al. 

1997) , B3 (Ficarra et al. 1985), and NVSS (Condon et al. 

1998) . Where the available radio maps make it possible, the 
flux densities of the lobes used to calculate their luminosi- 
ties arc cleared of the hot spots' emission. Columns (6) and 
(10) give the derived luminosities calculated using Ho=71 
kms"^Mpc"S f2m=0.27, and f2vac=0.73. 



4 DYNAMICAL AGE AND OTHER PHYSICAL 
PARAMETERS 

4.1 Independent solutions for the individual lobes 

In the first step of modelling the age and other physical 
parameters of the investigated GRGs, we fit them indepen- 
dently for each of the two lobes of a given source using their 
observational data given in Table 2 and the values of other 
free parameters of the model from Table 1, in paxticular with 
13=1.5. The dependence of the resulting model solutions of t, 
Qjet, po, and Qjet x t on ainj for the N-lobes of J0912-I-3510 
and J1155-I-4029, i.e. for these lobes for which we assumed 
6 < 90°, is shown in Fig. 2. The vertical line and the arrow 
indicate the value of CKinj that corresponds to the minimum 
of the kinetic energy delivered to the given lobe by the jet. 
The parameter values resulting from the fits are listed in 
columns (3) (7) in Table 3. 

As might be expected, the DYNAGE solutions of age 
and other physical parameters for the opposite lobes of the 
same source appear somewhat different. The question is 
whether these differences are statistically significant. To an- 
swer this question, we estimate an error on the fitted values 
of a given parameter. The error estimates for the fitted val- 
ues of ainj and age (t) of the lobes are included in Table 3. 
These errors indicate that a difference between the ages of 
the opposite lobes is always insignificant; however there is a 
marginal tendency of the larger lobes to be older than the 
opposite shorter ones. It is also worth noting that in at least 
three of the ten sources, a difference between the values of 
the effective initial spectral index Qinj is larger than 2.5-3 a. 
Such differences are possible if an evolution of the magnetic 
field, various energy losses and acceleration processes of the 
relativistic particles, and mixing of plasma at diflFerent ages 
- arc different at the heads of the opposite lobes. 

The above suggests that such differences in the ageing 
properties of the opposite lobes, especially in giant-sized ra- 
dio sources, may be caused by different environmental con- 
ditions. This is likely confirmed by the evident differences in 
the fitted values of po. Although this parameter denotes a 
density of the central radio core, in the DYNAGE algorithm 
we fit the term poa^ which appears in Eqs. (1), (2), and (3). 
If the values of oo and 8 are fixed, the fitting procedure 
gives the values of po slightly or significantly different for 
the opposite lobes. However, we can expect a single value 
of Po (and Qjet) for a given source. Therefore, in the next 
subsection we describe and check a 'self-consistent' solution 
for the investigated lobes. 
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Figure 1. a) Qjct — PO diagram for the NW lobe of J1313+6937. Tlic radio luminosities at the four frequencies indicated and used to 
fit the model are given in Table 2. The a^^j value for which the jet kinetic energy achieves a minimum is given, b) 'Goodness' of fit 
quantified by the A measure vs. the age for the same value of ainj . Its minimum indicates the age solution for the given lobe, of. Table 3. 



4.2 Self-consistent solution for the opposite lobes 

In this second step of the modeHing we averaged the values 
of Qjct and po found from fits for the opposite lobes (given in 
columns (5) and (6) of Table 3), and now treat them as the 
fixed free parameters of the model, (Qjet) and (po), respec- 
tively. Given these values, we can determine a value of /? for 
each of the two opposite lobes, hereafter denoted as /3s. c.. In 
order to do that, we assume that a power- law density profile 
of the ambient environment can be extended to distances as 
large as lengths of the lobes of GRGs. Although such an as- 
sumption can be invalid for distances larger than a few hun- 
dred of kpc (cf. Gopal-Krishna & Wiita 1987; Machalski et 
al. 2008), the ambient densities, pa, calculated for the inves- 
tigated GRGs from equation (1) with r = D are still within 
an acceptable range of density (~ 3 x 10~^® , ~ 2 x 10~^® 
kgm~^). Thus, equalising pa values in the independent solu- 
tion and the self-consistent solution, we have from equation 
(1) 

o ^ Iog((po)/Pa) 

and transforming equation (2) we calculate an expected age 
of a given lobe (in the frame of the self-consistent solution) 

from 

VcJ \^ (Qjet) J 

Finally, seeking a preservation of the minimum energy 
conditions, we search for the relevant value of ainj , hereafter 

denoted as alnj'^''. The resulting values of fls.c, Qinj'^ ^ ^s.c, 
and (tih/c)s.c. which is an average expansion speed of the 
lobe derived from the self-consistent solution, are given in 
columns (8)-(ll) of Table 3. 



5 DISCUSSION OF THE RESULTS AND 
CONCLUSIONS 

5.1 Ages and expansion speeds of the lobes 

The DYNAGE fits imply that (i) the formal ages of the op- 
posite lobes are somewhat different, though the difference 
is not large, given the errors of the fit; (ii) there is a weak 
trend in the sense the larger lobe is 'older', and (iii) the av- 
erage expansion speed of the larger lobe is usually higher 
than that of the shorter one. The implication (ii) is consis- 
tent with predictions of the simplified kinematic model for 
the jets' propagation (e.g. Longair & Riley 1979): if heads 
of the jets move through a uniform environment at the same 
speed at an angle 6 to the line of sight, the shorter lobe 
will appear younger by the time delay of (Di -|-D2)cot(6')/c. 
For the projected linear size Di + D2 = 1 Mpc and 6 ss 70°, 
an age difference will be less than 1.2 Myr. Therefore the 
derived age differences for the opposite lobes, much higher 
than 1-2 Myr, cannot be related to the kinematic effects 
only, but they likely indicate actual different jet's prop- 
agation conditions in the opposite directions through the 
galactic and/or intergalactic medium. The implication (iii) 
will be expected if both the lobes arc of the same age. In- 
deed, in spite of the formal age differences, in seven out of 
ten GRGs a higher expansion speed is found for the larger 
lobe. A faster expansion, in turn, may imply a thinner en- 
vironment and/or a higher jet power. Surprisingly, this is 
not the case for the four sources: J0912-I-3510, J1155-I-4029, 
J1313-F6937, and J1702-I-4217; the core density fitted for 
their larger lobes is higher than that for the shorter ones. 
Only the shorter (and formally older) lobes of the sources: 
J0927-I-3510, J1343-I-3758, and J1604-f3438 seem to expand 
slower due to a denser ambient medium at their side, how- 
ever the speed differences are insignificant, being of an order 
of the speed errors of about 0.005c-0.008c. 

The self-consistent solutions do not change the above 
picture radically. A median of the age quotient (older to 
younger) of about 1.2 is similar in both the independent so- 
lutions and in the self-consistent solutions. Still the expan- 
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Table 2. Observational parameters of the sources and their lobes 



lAU name 


z 


Larger 


D(kpc) 


u 


logP^ 


Shorter 


D(kpc) 


V 


logP^ 


Other name 




lobe 


R,x 


(MHz) 


(WHz-isr-i) 


lobe 


R,x 


(MHz) 


(WHz-lsr-l) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J0912+3510 


0.2489 


S 


901±19 


7C151 


24.881 


Nf 


578±82 


7C151 


24.556 








5.2±1.5 


G606 


24.436 




4.1±1.1 


G606 


24.154 










V1400 


24.156 






V1400 


23.892 










V4860 


23.710 






V4860 


23.456 


J0927+3510 


(0.55) 


SE 


1119±20 


7C151 


25.395 


NW 


1087±32 


7C151 


25.253 








4.6±1.2 


G606 


25.005 




4.4±1.2 


G606 


24.788 










V1400 


24.711 






V1400 


24.486 










V4860 


24.236 






V4860 


24.000 


Jl 155+4029 


(0.53) 


sw 


956±12 


G241 


25.346 


NEJ 


628±10 


7G151 


26.343 








8.1±1.4 


G605 


24.968 




3.2±0.8 


B408 


25.925 










V1400 


24.603 






V1400 


25.383 










V4860 


24.036 






V4860 


24.821 


J1313+6937 


0.1064 


SE 


422±30 


7C151 


24.806 


NW 


323±29 


7C151 


24.942 


DA340 






3.3±1.1 


G605 


24.389 




2.4±0.8 


G605 


24.516 










V1425 


24.095 






V1425 


24.229 










V4873 


23.593 






V4873 


23.759 


J1343+3758 


0.2267 


NE 


1404±36 


7C151 


24.409 


SW 


1059±11 


7C151 


24.716 








5.4±1.3 


W325 


24.200 




3.3±0.7 


W325 


24.423 










V1400 


23.811 






V1400 


23.981 










V4860 


23.290 






V4860 


23.382 


J1453+3308 


0.249 


S 


756±12 


G240 


24.953 


N 


570±12 


G240 


25.170 


B1450+333 






4.8±1.6 


G605 


24.704 




2.8±1.0 


G605 


24.933 










V1400 


24.356 






V1400 


24.586 










V4860 


23.705 






V4860 


24.029 


J1604+3438 


0.2817 


E 


443±17 


G239 


24.692 


W 


403±21 


G239 


24.638 








2.6±0. 


G614 


24.361 




2.4±0.4 


G614 


24.300 










G1265 


24.198 






G1265 


24.107 










V4860 


23.603 






V4860 


23.526 


J1604+3731 


0.814 


NW 


679±23 


G334 


25.708 


SE 


667±15 


G334 


25.804 


7C1602+376 






2.4±0.5 


G613 


25.439 




2.8±0.5 


G613 


25.530 










G1289 


25.219 






G1289 


25.318 










V4860 


24.544 






V4860 


24.697 


J1702+4217 


0.476 


SW 


663±47 


W325 


25.213 


NE 


497±47 


W325 


25.364 


7C1701+423 






3.2±0.7 


G602 


25.002 




1.9±0.5 


G602 


25.144 










V1425 


24.685 






V1425 


24.827 










V4860 


24.131 






V4860 


24.313 


J2312+1845 


0.427 


NE 


544±34 


V74 


26.751 


SW 


512±28 


V74 


26.843 


3C457 






4.1±1.2 


G334 


26.177 




3.4±1.1 


G334 


26.284 










V1425 


25.615 






V1425 


25.738 










V4866 


25.087 






V4866 


25.213 



Notes. A letter preceding the observing frequency indicates the radio telescope (array) used: G - GMRT, V - VLA, W - WSRT, or the 
surveys: 7C, B3 (cf. the text), f indicates the lobe with an inclination angle of 70° assumed; J indicates the lobe with an inclination 
angle of 50° assumed. The corresponding values of Re give the ratio of deprojected length of the lobes. 



sion speeds tend to be higher for the larger lobes. Again the 

exception seems to bo J1604+3731. The values of /3s. c. found 
by the fits are within an acceptable range (1.38, 1.76). It is 
worth noting that these values are always less than 2 which 
is the necessary condition for forming the head of the jet and 
to observe a source of the FRII type. On the other hand, the 
self-consistent solutions require a larger difference between 

(s c ) 

the effective injection spectral index ajj^j for the opposite 
lobes, than those found in the frame of the independent 
solutions. The question whether it really reflects different 
physical conditions governing the initial energy distribution 
of the rclativistic particles at the head of the lobes, or it 
is mostly related to an uncertainty of the fit and/or wrong 
assumptions in the model, is open. 

Another alternative self-consistent solution is plausible 
in which the age of the opposite lobes has to be the same. 



and any differences between their size and luminosity are 

due to an inhomogcncity (asymmetry) in density distribu- 
tion of the surrounding gaseous environment. Again, there is 
rather no physical circumstances for significant differences 
between the jet power and the central core density in the 
opposite directions along the jets' axis. In such a scenario, 
cither the same values of po, Qjct, and t would be assumed 
for both lobes or only postulated, while the values of ao, /3, 
and ainj, different for the two lobes, are determined by the 
fit. Rearrangement of equation (2) and substitution of D for 
L{t) gives 

Applying this equation separately for either lobe of a source 
and demanding equality of their age, jet power, core density. 
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Figure 2. The dependence of jet power, Qjet, age of the lobes, t, central core density, po, and kinetic energy delivered to the lobes 
during the time t, Qjet x t on the effective initial spectral index, ainj, for the N- lobes of two different GRGs from Table 2. The minima 
of the kinetic energy corresponding to the 'best solution' of the age and other physical parameters of the given lobe are marked with the 
arrows. The vertical dotted lines indicate proper y-intercepts of each of the four functions y = /(ainj)- 



Table 3. Fitted physical parameters of the lobes. Columns from (3) to (7) give results of the independent solution, while columns from 
(8) to (11) give results of the self-consistent solution 



Source 


Lobe 


Q^inj 


t 


log Qjet 


logpo 


Vh/c 


A.c. 




ts.c. 


(Vh/c)s.c. 








(Myr) 


(W) 


(kgm-3) 






(Myr) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


J0912+3510 


S 


0.493±.006 


90+17 


37.80 


-23.30 


0.033 


1.45 


<0.4 


(98) 


(0.031) 




N 


0.487±.007 


75+14 


37.47 


-23.51 


0.025 


1.57 


0.590 


66 


0.031 


J0927+3510 


SE 


0.490±.005 


41+4 


.38.55 


-23.99 


0.089 


1.52 


<0.4 


(44) 


(0.083) 




NW 


0.498±.005 


50+5 


38.38 


-23.90 


0.071 


1.18 


0.550 


46 


0.077 


J1155+4029 


sw 


0.570±.030 


130+35 


38.23 


-22.00 


0.024 


1.47 


0.632 


111 


0.028 




NE 


0.615+.019 


118+28 


38.55 


-22.13 


0.017 


1.54 


0.569 


133 


0.015 


J1313+6937 


SE 


0.508+.011 


152+22 


37.18 


-22.66 


0.010 


1.46 


0.486 


153 


0.010 




NW 


0.511+.011 


131+18 


37.15 


-22.80 


0.008 


1.55 


0.523 


131 


0.008 


J 1343+3758 


NE 


0.461±.018 


64+11 


37.90 


-24.24 


0.071 


1.58 


<0.4 


(69) 


(0.066) 




SW 


0.516±.012 


94+13 


37.76 


-23.94 


0.036 


1.44 


0.547 


87 


0.039 


J 1453+3308 


s 


0.563it.006 


173+23 


37.70 


-22.35 


0.014 


1.40 


0.550 


156 


0.016 




N 


0.518±.013 


127+17 


37.69 


-22.92 


0.015 


1.71 


0.544 


136 


0.014 


J 1604+3438 


E 


0.506±.007 


86+12 


37.40 


-23.40 


0.017 


1.50 


0.440 


90 


0.016 




W 


0.511+.008 


93+9 


37.29 


-23.39 


0.014 


1.50 


0.547 


89 


0.016 


J1604+3731 


NW 


0.571+.008 


54+8 


38.62 


-23.48 


0.042 


1.45 


0.590 


51 


0.045 




SE 


0.536+.009 


41+8 


38.73 


-23.70 


0.054 


1.57 


0.512 


43 


0.051 


J1702+4217 


SW 


0.568+.012 


74+15 


38.06 


-23.07 


0.030 


1.38 


0.544 


85 


0.026 




NE 


0.538+.012 


62+10 


38.11 


-23.73 


0.027 


1.76 


0.544 


70 


0.024 


J2312+1845 


NE 


0.567±.012 


86+21 


38.58 


-22.10 


0.022 


1.47 


0.574 


83 


0.023 




SW 


0.567±.013 


76+14 


38.64 


-22.23 


0.022 


1.54 


0.561 


79 


0.021 



and core radius, we have another equation involving two 
unknown quantities: /3i, and /32, where the fraction po/<3jet 
is eliminated 

„(/3l-fe) ^ (g2/Cl,2)<^-^^> 
° (Dl/Cl,l)(5-/31)- 

The values of two latter parameters depend on the energy 
distributions of particles injected into the opposite lobes and 
described by ainj parameters and on their observed luminosi- 
ties. This dependence cannot be expressed analytically, as 
explained in Section 2. Besides, the values of po and Qjet 
(as well as ainj ) cannot be determined independently of the 
unknown values of ao and (3. Therefore, such an 'alternative 



self-consistent solution' obviously cannot give any explicit 

results. 

In order to prove this we attempted to find such a solu- 
tion for the lobes of J0912+3510, the sample source with a 
rare asymmetry where the larger lobe is also much more lu- 
minous than the shorter one. Assuming for both lobes t = 82 
Myr, which is the average of the values determined in the 
'independent solution' (cf. Table 3), we fit the values of Qjet 
and Po for a number of combinations of the model param- 
eters ao, and Oinj. The results are shown in Fig. 3, where 
the abscissa axis gives the core density po scaled to the core 
radius ao=10 kpc. Exploring the parameter space we notice 
that (i) none of physically acceptable values of these param- 
eters can provide comparable values of Qjet for the opposite 
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lobes (cf. the similar paper of Brocksopp et al. 2007 and 
their Fig. 4), and (ii) allowing different values of Qjet and 
po for either lobe, the assumed age value can be achieved 
with any combination of ao, fi and ainj parameters. The 
above calculations confirm that there is no unique 'alterna- 
tive self-consistent solution'. 

There are two possible causes of different 'dynamical 
age' solutions of the opposite lobes of a FRII-type source, 
especially a GRG in which the effect of projection is likely to 
bo neglected. The first cause, the most probable one, could 
be related to an unknown actual distribution of magnetic 
fields, which may depart from the equipartition conditions 
in a different way in cither lobe. This is confirmed by X-ray 
imaging observations which allow measurements of electron 
energies of radio lobes and magnetic fields. The derived elec- 
tron energy densities often exceed those of magnetic fields by 
a factor of a few to several dozen. For example, the ratio ^ of 
about 0.16 has been found in the lobes of the galaxy CenB 
(Tashiro et al. 1998) and of 0.20 and 0.13 in the E lobe and 
W lobe of the galaxy For A (NGC 1316), respectively (Isobo 
et al. 2006; Tashiro et al. 2008). Our calculations show that 
increasing or decreasing the magnetic-field energy density, 
Mb, in respect to the equipartition condition (given in equa- 
tion (4)) results in an increase or decrease of the fitted age, 
respectively, which agrees with the finding of Parma et al. 
(1999). For example, the four- fold increase of the ratio f in 
the N-lobe and the corresponding decrease in the S-lobe of 
J0912+3510 results in about 10%-15% increase or decrease 
of their fitted age, respectively, i.e. eliminates the age dif- 
ference found with the 'independent solution'. This is worth 
to emphasize that this variation of ^ is not equivalent of 
changing ainj. In both lobes a minimum of Qjet x t is still 
around 0.49, which strongly supports the conclusion that the 
minimum of the jets' kinetic energy appearing in the age so- 
lutions is not an artefact of numerical calculations, but has a 
real physical meaning. The corresponding effective spectral 
index ainj is always within the narrow range of p ~ 2.0-2.4 
suggested by the 'non-relativistic shock' paradigm for the 
origin of non-thermal electrons within the heads of FRII- 
type radio sources (e.g. Blandford & Eichler 1987; Heavens 
& Meisenheimer 1987). 

The second cause, less probable one, could be related 
to a multi-episod jet activity, strongly suggested by observa- 
tions of double-double radio galaxies (DDRGs, e.g. Kaiser et 
al. 2000; SaripalU et al. 2002; Brocksopp et al. 2007). If the 
jet activity has been interrupted and its multiple episodes 
happened in opposite directions by chance, an 'effective' age 
of the lobes would be different. There is, at least, one piece of 
observational evidence for such a behaviour, i.e. the northern 
middle lobe of the galaxy CenA which has no counterpart 
in its southern lobe (Morganti et al. 1999). Its spectral age 
of about 20-30 Myr, estimated by Hardcastle et al. (2008), 
is much higher than the time scale of about 10* yr for in- 
stabilities in the accretion processes. 

5.2 A relation between the spectral age and the 
dynamical age 

The spectral age distribution along the jets' axes of the ten 
investigated GRGs was analysed in Paper II. In this subsec- 
tion we compare (i) the initial spectral indices resulting from 
the DYNAGE and SYNAGE fits, and (ii) the dynamical ages 



Of 



0.60 




0.56 




0.52 




0.48 








0.60 , 


N-lobe 










48 " 













-23.8 -2'A.l -2'.i.0 -22.6 

log(poja„[kpc]/10!'') 

Figure 3. Model solutions for the lobes of J0912+3510 at the 
age of 82 Myr for both lobes. Abscissa gives the (log) core density 
scaled to the core radius of 10 kpc. The areas confined by dotted 
lines indicate very large parameter space within which both lobes 
would have the assumed age. 

derived in this paper and the synchrotron ages analysed in 
Paper II and supplemented in this paper. As described in 
the Introduction, the notion of a spectral age of the source 
(or its lobes) is not explicit; it is usually considered as an 
age of the emitting particles which can differ in different 
parts of the source. Therefore in Paper II we derived two 
estimates of this age: the age of the oldest detected particles 
in the regions of the observed emission from the lobes, and 
the age values resulting from extrapolation of a linear regres- 
sion of the age on the distance (i.e. a characteristic speed, 
Vsep, which is an indication of the speed of the lobe material 
relative to the hot spots) to the radio core. Besides, both 
these estimates were derived using two different equiparti- 
tion magnetic field formulae, the classical formula of Miley 
(1980) and the revised formula of Beck & Krause (2005). In 
each case the spectral age was computed using a spectral 
break frequency, Vhr, resulting from the SYNAGE fit of the 
Jaffe & Perola (1973; JP) model to the radio spectrum of a 
given lobe for which the JP fit was mostly the best compared 
with those of the Kardsishev-Pacholczyk (Pacholczyk 1970; 
KP) or 'continuous injection' (Kardashev 1962; CI) models. 

The DYNAGE method for estimating the dynamical 
age of FRII-type radio sources is based on the analytical 
KDA model which assumes a continuous delivery of energy 
to the lobes through the jets with a constant power. There- 
fore, in order to compare the initial spectral indices resulting 
from the DYNAGE and SYNAGE fits, we should take into 
account the SYNAGE indices fitted with the CI model. How- 
ever, this model of energy losses does not account for an adi- 
abatic expansion of the cocoon, while KDA and DYNAGE 
do. Therefore, a model much closer to DYNAGE is the CIE 
model of Murgia (1996) which accounts for an adiabatic ex- 
pansion of the volume of source (lobe, cocoon) and a related 
evolution of the magnetic field strength in the form r{t) oc 
{V{t) oc **'=) and B{t) oc t"". Though each SYNAGE fit of 
the CIE model to the data points in the radio spectrum of 
a given lobe was the worst as compared with those of the 
JP, KP, and CI models - the resulting spectral ages, tcie, 
should be be closer to the dynamical ages estimated in this 
paper. 
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The break frequency for the CIE model is i^br.ciB = 
(fe - 2m - l)^i4>r,ci- In the SYNAGE algorithm, where k=l 
and m=— 2, i^r,ciE = 16 iA)r,ci are expected for a spectrum 
of same initial energy distribution, the same age, and fi- 
nal magnetic field. However, for 20 lobes of the analysed 
GRGs the ratio i^r.ciE/i^br.ci is consistently much less than 
16, in the four sample sources is even less than 1. This is 
caused by a frequent large difference between aci and aciE 
fitted as a free parameter in SYNAGE. The values of low- 
frequency spectral index and break frequency for the in- 
vestigated lobes, fitted with SYNAGE in the frame of CI 
and CIE models, as well as the values of corresponding syn- 
chrotron ages are given in Table 4. The value of magnetic 
field strength Bdyn, derived from the magnetic energy den- 
sity in this paper (cf. equations (4) and (5)) and used to 
calculate the ages, is given in column (3) . The values of aci , 
I'br.ci, aciE, and iA,T,cm are given in columns (4), (5), (7), 
and (8), respectively. The relevant ages, rci and tcie, given 
in columns (6) and (9), respectively, arc calculated from 

^ci = 50.3-g5^J^{i.b,,ci(l + z)}-'^' [Myr], (6) 
and 

TCIE = 201.2 „f {iA.r,ciE(l + [Myr], (7) 

where Bic[nT]=0. 318(1 + z)'^ is the inverse- Compton mag- 
netic field strength equivalent to the cosmic microwave back- 
ground radiation. Finally, the differences between the fitted 
low-frequency slope and the injected spectral index deter- 
mined with DYNAGE, aci — ounj , and qcie — ainj arc given 
in columns (10) and (11), while the ratios between the dy- 
namical age t derive in this paper and the synchrotron ages 
TCI and Tcie - in columns (12) and (13), respectively. 

The data in Table 4 show that the equipartition mag- 
netic fields derived with the DYNAGE algorithm are usually 
stronger than those calculated with the classical Miley's for- 
mula. On the second hand, all these fields are weaker that 
those calculated with the revised formula of Beck & Krause 
(2005) published in Paper I. The reason for the latter efi^ect 
is that the observed curved radio spectra were integrated, 
while Beck & Krause assumed that the proton spectrum is 
straight and any steepening in the radio spectrum is due to 
energy losses of the electrons. If, however, the proton spec- 
trum is similarly curved as the electron spectrum, the Beck 
& Krause values would be too high. Nevertheless, the result- 
ing spectral ages are not too much sensitive to the magnetic 
field strengths, depending more on the spectral break i\,t 
and on the ratio B/Bic- The examples of the spectral best 
fit for the lobes of two different sample sources are shown 
in Fig. 4. Three different fit results for the N lobes of the 
sources J0912-I-3510 and J1155-I-4029 are compared. Note 
that the quality of the fits, expressed by the reduced val- 
ues, are comparable within the frequency range for which 
fiux-density data have been available. The fits clearly show 
that the frequency range crucial for a better discrimination 
between the models is below 100 MHz, i.e. the range where 
near-future observations with the LOFAR array will be de- 
cisive. 

The data in Table 4 also exhibit unexpected discrep- 
ances between the fitted values of aci (or acm) and i^br 
in opposite lobes of some sample sources, e.g. J1155-I-4029 



J1155+4029 




u [MHz] 



Figure 5. Radio spectra of the opposite lobes of J1155-|-4029 and 
their spectral curvature parameters fitted with the original CIE 
model. 



and J1343-I-3758. The most discrepant are the values of 
t'br.ciE in J1155-I-4029. This is shown in Fig. 5 where the 
observed spectra of its opposite lobes are compared. Both 
slopes and spectral shapes look very similar, nevertheless 
the formal fitted values of ainj and especially of ;4,r,ciE are 
dramatically different. How much these values are uncer- 
tain is indicated by their errors provided by SYNAGE: 32 
MHz< i^br.ciE <241 GHz and 0.08 MHz< J^br.ciE <1-18 
GHz for the SW lobe and NE lobe, respectively. This is 
worth to notice that J1155-I-4029 has the largest asymme- 
try in the lobes' luminosities. 

The fit results indicate that a difference between the 
dynamical age, t, and the synchrotron age, rci or tcie is 
related to a difference between the values of ainj (column 
(3) of Table 3), oci or aciE (columns 4 and 7 of Table 4). 
As this was noted in paper of Machalski et al. (2007), 
Qinj values derived with DYNAGE are usually lower than 
those fitted with SYNAGE. An opposite situation is rare; 
an example is the NE lobe of J1155-I-4029 for which the 
SYNAGE fit gave qci=0.537±0.021 and qcie=0.5801o:o5 
while Qinj=0.615±0.019 is found in this paper. Such a case 
is likely caused by a very uncertain observational data at 
low frequencies. In general, an increase of the difference 
Usyn — Qinj corresponds to an increase of the ratio t/rsyn, 
where the label 'syn' indicates either aci or aciE. This cor- 
relation (shown in Fig. 6) is strong; the correlation coeffi- 
cient between the above two variables for our twenty lobes 
is -1-0.837 in the case of t/rci and -1-0.884 in the case of 
t/Tcm- The data in Table 4 show that in 13 of 20 investi- 
gated lobes the synchrotron age calculated with fc=l and 
m=— 2 is higher than the dynamical age determined in this 
paper (the ratio t/rcm is less than unity in Fig. 6). Indeed, 
the corresponding factors in the KDA model of the jet's 
dynamics, given by equations (2) and (5), are fc=6/7 and 
m=— 11/14 for the assumed values of /3, Fb, and Fc. As a 
result, we rather would expect tcie = (10/7) tci instead of 
Tcie = 4 tci- In fact, in Fig. 6 the y- intercepts of t/rci and 
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Table 4. Dynamical and synchrotron ages of the lobes 



Source 
(1) 


Lobe 
(2) 


(nT) 
(3) 


"CI 
(4) 


(GHz) 
(5) 


TCI 

(Myr) 
(6) 


a 


CIE 

(7) 


t'br.CIE 

(GHz) 
(8) 


TCIE 
(Myr) 

(9) 


"CI- 
"inj 

(10) 


"CIE- 
Oinj 
(11) 


t/ 

TCI 

(12) 


t/ 

TCIE 

(13) 


J0912+3510 


S 





16 





717 


18 


87 


15.3 





528 


4.97 


119. 





224 





010 


5 


88 





76 




N 





16 





637 


8 


21 


23.1 





503 


4.91 


120. 





150 





016 


3 


24 





63 


J0927+3510 


SE 





18 





684 


3 


74 


14.4 





541 


1.63 


87.2 





194 





051 


2 


85 





47 




NW 





14 





434 





30 


45.7 





554 


1.77 


75.3 


-0 


064 





056 


1 


09 





66 


J1155+4029 


SW 





34 





914 


15 


36 


15.3 





758 


4.29 


68.4 





344 





188 


8 


47 


1 


90 




NE 





62 





537 





11 


103. 





580 


(0.046) 


(636) 


-0 


078 


-0 


035 


1 


15 





19 


J1313+6937 


SE 





23 





613 


2 


54 


70.4 





490 


0.82 


496. 





105 


-0 


018 


2 


16 





31 




NW 





28 





643 


4 


29 


116. 





477 


0.72 


519. 





132 


-0 


034 


1 


13 





25 


J1343+3758 


NE 





08 





571 


2 


31 


35.9 





502 


1.59 
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Figure 4. Exemplary radio spectra of the lobes and fits of the three models of energy losses: CI (a pure continuous injection) — the solid 
line, original CIE (continuous injection and expansion - the dashed line, and DYNAGE the dotted line. The values of ainj fitted for 
each of the models and the corresponding reduced values are given for a comparison. 



t/rcm at Oayn — ainj=0 differ by ~ 2, not by 4. Therefore, we 
argue that the CI model underestimates the lobe's age sinee 
it does not take into account the expansion losses, while the 
CIE one overestimates the lobe's age since the original ex- 
pansion factors (fc=l and m=— 2) are likely too extreme. In 
between of these, the DYNAGE age solution should be bet- 
ter since the expansion parameters are connected to actual 
geometry of the lobes for each specific GRGs. 

The plot in Fig. 6 shows that the dynamical age of the 
GRG's lobes is about 1 to 5 times larger than their raxliar 
tive age. Very similar age ratios were found by Parma et 



al. (1999) for much smaller, low luminosity radio galaxies. 

Beside the causes of this effect already discussed in the liter- 
ature, there are two other possible explanation of this effect: 

(1) As expected, both the synchrotron and the dynam- 
ical ages are sensitive to the o^inj parameter. However, it 
seems that an ainj value determined with the DYNAGE 
method from the minimum kinetic energy of the jets is less 
sensitive to an uncertainty of the low-frequency spectrum of 
a source than values of aci or aciE found with SYNAGE. 
The only difference between the spectra provided with the 
model accounting for the adiabatic losses and the model not 
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Figure 6. Ratio of the dynamical age of the lobes, t, and their 

synchrotron age, rsyn, vs. the difference between asyn — ainj where 
Qsyn means either aci or i^ciE- The circles indicate the ratio of 
t and TQi while the crosses - the ratio of t and tcie- The two 
dotted lines show linear regression lines of the above ratios on 
the asyn — cuinj coordinate. 



accounting for them - is a larger range of transition of the 

spectral slope from Qinj to Oinj + 0.5 in the former model. 
Therefore, in further discussion given below, we do not dif- 
ferentiate between aci and aciB. The ainj and aci parame- 
ters both determine the slope of the spectrum at frequencies 
where radiative losses are not yet important. 

In the classical spectral ageing analysis the CI model in- 
cludes the continuum injection of relativistic particles having 
all the time the same power-law energy distribution. In the 
presence of energy losses due to magnetic fields the observed 
spectrum, being a sum of the emission from the various par- 
ticle populations at different synchrotron ages, is charac- 
terized by a curvature between its low-frequency and high- 
frequency ends. This curvature is formed by many spectra 
with different normalizations and break frequencies. Implicit 
in this model is the assumption that the field strength is the 
same in the region of injection as in the outflow region. But 
the DYNAGE model is different because it accounts for the 
evolution of the magnetic field strength in a different way 
than the original CIE model does that. Hence the spectra 
calculated with it have different curvature. Probably the ac- 
tual spectra have some curvature in them as well. SYNAGE 
cannot cope with such a curvature, thus fits a spectrum, 
single-broken between values of an aci and aci + 0.5, to the 
data points. However, because the observed spectrum is al- 
ready affected by radiative losses below the fitted break, aci 
(and qcie) comes out rather steep. Contrary to that, DY- 
NAGE can cope with the above curvature because, basing on 
a more realistic model of dynamical expansion of the source 
(cocoon), it can better predict its spectrum. Therefore it can 
fit an aged spectrum with a fiatter ainj , because the part of 
the spectrum already affected by radiative losses extends to 
much lower frequencies than in the case of the classical ap- 



proach. However, it should be noted that the KDA model for 
the cocoon's emission is based on the approximation of the 
synchrotron kernel radiation with a 'delta'-function. This 
approximation could artificially enhanced the spectral cur- 
vature, in the sense that it could result much more pro- 
nounced than it really is. 

Resuming, we argue that usually DYNAGE, finding a 
flatter ainj, is also less sensitive to its exact value. In other 
words, DYNAGE can afford radiative effects at lower fre- 
quencies than SYNAGE. Radiative effects at lower frocjuen- 
cies means older age, thus one can suspect that this effect 
is stronger at older ages than at younger ages where the 
spectra tend to be 'straightcr', i.e. with only a single break 
and less curvature. But our, though very limited, statistics 
shows something else. The difference aciE — ainj seems to 
depend on redshift but not on the age. The correlation co- 
efficient in the correlation between this difference and the 
redshift is -1-0.63. Again, this effect can support our thesis 
that DYNAGE can better afford radiative losses at low fre- 
quencies, where an intrinsic spectral curvature is shifted by 
the redshift, than the classical spectral ageing analysis. 

(2) There is another factor having an influence on the 
difference between t and r. This is the axial ratio Rt- 
Note that in the classical energy equipartition formula, 
Beq oc {L/Vcf^'', the cocoon's (lobe's) volume is completely 
insensitive of a value of Rt what originate from the spher- 
ical geometry assumed in the classical formula of Pachol- 
czyk (1970). Oppositely, because of the cylindrical geome- 
try assumed in KDA and DYNAGE, the cocoon's volume is 
strictly dependent of Rt which is K oc R^^ (its length D 
is measured much more precisely than the transversal base 
diameter). A variation of Rt results in a change of the pres- 
sure ratio Vhc which determines the lobe (cocoon) energy 
density. This, in turn, specifies the magnetic field strength 
via the magnetic energy density (cf. equations (3), (4), and 
(5)). The data in Tables 3 and 4 confirm that the smaller 
the differences aciE — ainj (and closer to unity the ratio be- 
tween magnetic field strength calculated with the fornuila of 
Miley (1980) and the field determined with equation (5) in 
this paper, Bdyn), the closer the dynamical and radiative 
ages. 
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